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Abstract 


Background: In heavy ion collision at the Fermi energies Isospin equilibration processes occur¬ 
ring when nuclei with different charge/mass asymmetries interacts have been investigated to get 
information on the nucleon-nucleon Iso-vectorial effective interaction. 

Purpose: In this paper, for the system i8 Ca + 27 Al at 40 MeV/nucleon, we investigate on this 
process by means of an observable tightly linked to isospin equilibration processes and sensitive 
in exclusive way to the dynamical stage of the collision. From the comparison with dynamical 
model calculations we want also to obtain information on the Iso-vectorial effective microscopic 
interaction. 

Method: The average time derivative of the total dipole associated to the relative motion of 
all emitted charged particles and fragments has been determined from the measured charges and 
velocities by using the 47 r multi-detector CHIMERA. The average has been determined for semi¬ 
peripheral collisions and for different charges Z\, of the biggest produced fragment. Experimental 
evidences collected for the systems 2 ' Al+ 48 Ca and 2 ‘ Al+ 40 Ca at 40 MeV/nucleon used to support 
this novel method of investigation are also discussed. 

Results: The data analysis shows a clear signature of a trend to the global Isospin equilibration 
of the system for increasing differences of the Zj, values with respect to the projectile charge. 

Conclusions: The comparison with CoMD-II calculations gives the best agreement with data 
using a stiffness 7 value for the Iso-vectorial interaction in the range 7 ~ 1 — 1.2. Moreover, the same 
comparison allows to estimate the non-negligible contribution to the global isospin equilibration 
process given by the un-detected emitted neutrons. 

PACS numbers: 25.70.-z, 21.30.Fe 


* papa@ct.infn.it 
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I. INTRODUCTION 


Experimental evidences on Heavy Ion Collisions highlight, in different aspects, processes 
evolving on different time scales. At the Fermi energies semiclassical dynamical models can 
not describe the system during its overall time evolution. According to the displayed phe¬ 
nomenology, the collected data are usually described as produced by a fast pre-equilibrium 


stage described through dynamical models [1-3] and later stage processes described by statis¬ 
tical decay models { 4 ]. The statistical contribution is usually separated from the dynamical 
one through cut/extrapolation procedures applied to angular correlation and/or to particle 
kinetic energy spectra (see as an example [5j, 6|). When clearly identified 7], the attempt 
to measure observables in principle closely linked to only one of the two regime is therefore 
highly desirable. In fact, this eventuality allows to decouple effects related to the two classes 
of mechanisms that are linked to rather different nuclear matter macroscopic properties. In 
the last decades great efforts have been performed to extract information on the nuclear 
iso-vectorial forces by studying charge/mass asymmetric systems RR . These attempts 
concerns both the dynamical stage and the statistical decay of the produced hot sources. In 
particular in this last stage the isospin and excitation energy dependence of the level density 


formula play a key role and it is currently under investigation 


10, 12, 18 -[20]. 


A phenomenon closely linked with iso-vectorial-forces is the well known process leading 
to the redistribution in phase-space of the charge/mass excess /3 = 21] of the emitted 


particles and fragments (Z is the charge, N is the neutron number and A the mass num¬ 
ber). This phenomenon commonly referred as charge/mass or isospin equilibration process is 
rather complex in the Fermi energy domain, especially when finite size effects of the studied 
system have to be properly taken into account. The charge/mass distributions related to 
the final fragments and particles are affected by the pre-equilibrium stage, which includes 
particles and fragments production in the mid-rapidity region, prompt emission, transfer of 
mass and charge between the main fragments and finally particles/fragments emission from 
rot equilibrated sources through a multi-step statistical cascade. As an example, in Refs 


n],y, 


21-23] the transfer of mass-charge between the main partners has been investigated 


through the study of the isospin transport ratio obtained starting from the isotopic distribu¬ 
tions produced near the projectile rapidity for medium-heavy symmetric/quasy-symmetric 
systems. 
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In this paper we report on results of investigations on this equilibration process for the 
system 48 Ca + 27 Al at 40 MeV/nucleon starting from a different and/or complementary 
point of view. The measurement was performed with the CHIMERA multi-detector 24] at 
Laboratori Nazionali del Sud di Catania (Italy). The main goal of the experiment was to 
evaluate, for well reconstructed events, belonging to selected classes /C, the quantity: 


m 

0) = (J2 Z irfi-^c-m.))>C- 


( 1 ) 


2=1 


The brackets indicate the average value over the ensemble 1C. Zi,Vi,m are the charges, 
laboratory velocities, charged particle multiplicity respectively of the produced particles in 
the selected class of events, respectively. Finally T^ c . m . is the center of mass (c.m.) velocity. 
We note that in this expression the contribution of produced neutral particles is implicitly 
contained in V c .m.- The interest on this quantity was triggered by two main reasons: 

a) as shown in Refs 7, 25j this quantity is closely linked with charge/mass equilibration 
process because it represents the average time derivative of the total dipolar signal in the 
asymptotic stage (expressed in unit of e). In fact, as an example, for binary systems, in 
absence of dynamical neutron-proton collective motion we have: (l5) = = |(/r)((/3 2 ) — 

(A)X(^i - A)}. H is the reduced mass number of the system, /? i ,/ 3 2 are the isospin 
asymmetries of the two partners 1 and 2 and finally I* i and 2 the related velocities. In 
the above expression we have supposed negligible the correlation of fluctuations between 
charge/mass ratios of the partners and their relative velocity. On the other limit, the same 
quantity is zero if evaluated for a system represented by an equilibrated source before or 
after the statistical decay. As shown from dynamical microscopic calculations in a collision 
process between two nuclei having different charge/mass asymmetries, | (D) | changes during 
the time towards smaller values in the pre-equilibrium stage (spontaneous approach to the 
equilibrium) producing y-ray emission through the excitation of a more or less damped 


dipolar dynamical mode [7, 1261; 


b) because of the symmetries of the statistical decay mode, (i5) is not affected by the 
statistical emission of all the produced sources in later stages, as it is shown in Ref. ]7]. 
This essentially happens because, due to the vectorial kinematical character of this quantity 
, for well reconstructed events statistical effects are self-averaged to zero. Therefore (D) is 
a rather well suited global variable to selectively evidence dynamical effects related to the 
Isospin equilibration process. 
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II. THE EXPERIMENTAL PROCEDURE AND DATA ANALYSIS 


The experiment was carried out by using 40 MeV/nucleon 48 Ca and 27 Al beams, at the 
LNS Super-Conducting Cyclotron. The beam impinged on 400-^ 27 Al and about 1200-^ 
48 Ca, 40 Ca targets. The chosen combinations were the following: 48 Ca + 2 ' Al as the main 
system to be investigated, the 27 Al+ 48 Ca as the auxiliary system (see the following) and the 
charge/mass quasi-symmetric 27 Al-\- 40 Ca as reference system. Charges, masses, energies and 
velocities of the produced particles and fragments were measured with the An Multi-Detector 


CHIMERA 


24 


27| . In particular, the A E — E technique was employed for Z identification of 


fragments punching through the silicon detectors and additionally for isotopic identification 
of fragments with atomic numbers Z < 10. Mass identification is performed with the time- 
of-flight (TOF) technique by using the time signal from silicon detectors with respect to 
the time reference of the radio-frequency signal from the cyclotron. The TOF technique is 
basically used for velocity measurements of heavy ions. This technique is also essential for 
the mass and charge indirect determination of slow TLF’s (Target Like Fragments) stopped 
in the silicon detectors. Energetic light charged particles, stopped in the scintillator crystal, 
are identified by applying the ” fast-slow” discrimination method 28]. 


In the following we report results collected for events with a multiplicity of detected 
charged particles greater than or equal to 2. We have chosen rather restrictive selection 
criteria to identify the ’’good” reconstructed events. These conditions are suggested from 
calculations to obtain the ” invariance” of the investigated quantity with respect the statisti¬ 
cal decay mode (see Table I). For the main system we have selected events for which the total 
identified charge Zf ot = 33. Checks are being made to see at what extent the previous con¬ 
dition can be slightly released always keeping the above mentioned ’’invariance” condition 
to an acceptable level. The total detected mass was chosen in the interval 62 < Af ot < 78. 
The total measured momentum along the beam axis has been selected within 70% of the 
theoretical value (422 amu • cm/nsec). 


Analogous conditions have been imposed for the others investigated systems, taking into 
account the differences in the total mass and c.m. velocity. The well reconstructed events 
have been classified according to the charge of the biggest detected fragment Z\ and accord¬ 
ing to the estimated total kinetic energy loss TKEL = /j,E™ — YYLi \(MiV 2 — ^toYc.m.,d) ■ 
Where E™ is the incident energy per nucleon and /i is reduced mass number of the impinging 
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FIG. 1. (Color online) Charges Z of the detected fragments versus their velocities Vz along the 
beam axis are shown for the main system (panel a)) and the auxiliary one (panel b). 

0 2 4 6 8 10 


nuclei. In the above expression Mj,V) are the measured masses and velocities; Mf ot is the 
total measured mass and \^ c .m.,d the related c.m. velocity associated to the detected charged 
particles. Therefore, each event has been characterized also through these two last quan¬ 
tities, that in the TKEL evaluation can globally compensate the uncertainties due to the 
non perfect mass identification and velocity measurement. For the main system i8 Ca+ 27 Al 
and the auxiliary one 27 Al + 48 Ca, in Fig.l we show the charge of the detected fragments as 
function of their velocity Vz along the beam direction. In both cases the bi-dimensional plots 
show the dominance of processes producing TLF’s and PLF’s (Projectile Like fragments). 
The thin black curves represent the threshold of the A E — E charge identification technique 
as function of the fragments velocity evaluated for the typical Silicon detector thickness of 
the CHIMERA apparatus (about 300 fim). 


For each velocity value, fragments having a charge smaller than the value plotted through 
the curves can be directly assigned by means of the A E — E technique. Fragments having 
an higher charge will be stopped in the silicon detector an then the charge identification is 
obtained in an indirect way through the mass determination obtained by means of TOF and 


energy measurements and by using the Charity prescription 


29]. 
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FIG. 2. For TKEL<350 MeV, (panel a): Charge distribution of PLF’s punching through the 
Si detectors (direct estimate based on the A E — E method) for the 21 Al + 48 Ca system at 40 
MeV/nucleon (panel b): Charge distribution of TLF’s stopped in the Si detectors (indirect estimate 
based on the TOF method) for the 48 C*a + 2 ' Al system at 40 MeV/nucleon 

Therefore the comparison between the TLF charge distribution obtained in the main 
system in the indirect way (upper portion of the Z -Vz plot in Fig.l(a)) and the one related 
to the PLF for the auxiliary system, directly assigned, (lower portion of Z-Vz plot in Fig.l(b)) 
allows to evidence eventual systematic errors in the indirect charge assignment associated 
to the main system. The comparison has been performed for different windows of TKEL. 

It results that for TKEL lower than about 350 MeV the obtained charge distributions are 
similar and no systematic error°is evidenced. ThisAesult is shown in Fig.2 in panels a) and 
b) for the main system and the auxiliar one. For larger value of TKEL, the two distributions 
are different. In fact in this case due to the opposite kinematical conditions and to the finite 
geometrical efficiency of the detector, different reaction mechanisms are selected on average 
and the comparison is difficult. 

The present study is focused on the investigation of the behavior of the component of 
along the beam axis Z, Dz- As for the case of the TKEL determination and especially to 
compensate the eventual systematic errors in the fragment velocity determination, instead 
of the theoretical value we are necessarily induced to use in eq.(l) the value of 
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FIG. 3. (Color online) (panel a): For the system i8 Ca + 27 Al at 40 MeV/nucleon the measured 
values of D c z are plotted for different Zb associated to the selected events (charged multiplicity 
m > 2). The dot-dashed vertical lines indicates the reference limiting values D m (see the text), 
(panels b),c),d)): D c z distributions obtained as projection of the above bi-dimensional plot for 
different Zb intervals. In the panel b) the D c z spectrum for Zb =20 and for quasi-elastic events 
TKEL<70 MeV is plotted with star symbols. 

V r z m d evaluated, event by event, from the velocities along the beam axis of all the detected 
charged particles. The obtained quantity will be named in the following D z . D z therefore 
represents a partial dipolar signal related to the intrinsic motion of the subsystem formed 
by all the produced charged particle. The contribution related to the global relative motion 
of the undetected free neutrons is not included. In the next sections we will describe a way 
to estimate the average total signal associated to D z . In Fig.3 (panel a)) ,for the main 
system , we show the correlation plot Zb .v.s. D c z for the selected events. The ridge in the 
plot highlights an increasing trend of ( D z ) going from negative values to almost zero for 
decreasing values of Zb respect to Zplf • We note that according to the expression given for 
D m , in the initial configuration the system should exhibit a limiting value of ( D z ) (grazing 
collisions) close to about -9.8 cm/nsec. The increasing average values of ( D z ) for Zb different 
from projectile atomic number (less peripheral collision), represents a clear signature of the 
evolution through charge/mass equilibration values (see next section). The limiting value 
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corresponding to almost ’’grazing” collisions (along the beam axis), is evidenced in Fig.3 by 
dot-dashed vertical lines. Panels b),c) and d) show the projections of the bi-dimensional 
plot for different intervals of Z & and we can clearly see the trend of ( D c z ). In particular, 
in panel b) the spectrum with star symbols is obtained by imposing a value of TKEL<70 
MeV typical for quasi-elastic processes. The large fluctuations of D c z around the average 
value are due to physical reasons (the particular ’’history” of the each event) and to the 
measurement procedure reflecting the related uncertainties. In Fig.4 we show analogous 
plots for the reference isospin quasi-symmetric system 27 Al + 40 Ca at 40 MeV/nucleon. In 
this case the limiting value for ’’grazing” collision, D m is about -2.6 cm/nsec and, how can 
be clearly seen, the experimental plots shows values close to zero and an enhancement near 
the D m value (Fig.4b). The check on this system as compared to the main one ensure us on 
the good level of confidence obtained in the determination of ( D c z ). 


III. COMPARISON WITH COMD-II CALCULATIONS 


To obtain information on the behavior of the Iso-vectorial interaction starting from the 


present experimental data, we have performed CoMD-II 


30 


31] calculations in the interval 


of impact parameters 6 = 4 — 9 fm (the weight of each b is chosen proportionally to b 
itself). Dynamical calculations have been followed up to about 500fm/c. After this primary 
stage, the produced hot main sources have on average an excitation energy lower than 2.5 
MeV/nucelon. For each generated event, a second stage of a multi-step statistical decay 
has been simulated through the Monte Carlo GEMINI code jfkj . The obtained results have 
been analyzed with an implemented code that takes into account the main filtering effects of 
the experimental apparatus. This include geometrical acceptance of the identified particles 
and the main selection criteria used in the analysis of experimental data. 

In particular the CoMD-II calculations have been performed with stiffness parameter 
values concerning the Iso-vectorial interaction 7 = 0.5, 7 = 0.8, 7 = 1, 7 = 1.2, 7 = 1.5 
and symmetry energy about 32 MeV. According to previously performed investigations, the 
parameters of the effective Skyrme interaction corresponding to a compressibility of about 
220 MeV have been chosen following Ref. 33]. In Fig.5 (panel A),B) C)) we show the com¬ 
parison of the calculations with the measured value of (D z ) (red points) evaluated starting 
from the already shown correlation plot in Fig.3. The values are referred to different Z & 
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FIG. 4. (Color online) Same as the Fig .2 but for the reference system 27 Al + 40 Ca 

windows and TKEL<350 MeV. The vertical bars indicate the evaluated statistical errors 
which include also, conservatively, a non-systematic uncertainty on the measured fragment 
charge A Z = ±1. The corresponding theoretical values are plotted in the different panels 
with black symbols. The obtained values for 7 = 0.5 are rather different from the experi¬ 
mental results, they are out of the range value plotted in Fig.5. For this case a maximum 
value about 7 cm/nsec is obtained for Z h — 9 — 11 while the minimum value about -20 
cm/nsec is obtained for Z b = 21 — 22. The comparison of the others calculations with the 
experimental data is quantified in the figure through the shown y 2 value. According to 
this parameter, the case with 7 = 1 gives the best comparison even if it is only slightly 
better than the case 7 = 1.2. The behavior of (D c z ) (related to only charged fragments and 
particles) can be understood in a qualitative way according to the following considerations: 
as shown in the previous figure (apart from the case 7 = 0 . 8 ) for Z b smaller than Zppp 
(dominance of mass transfer from PLF to TLF) and for Z b larger than Zplp (on-set of 
fusion,incomplete fusion) ( D c z ) increases (there are on average more charge/mass symmetric 
produced fragments). In fact,in this case it can be verified that (V r z m d ) < V c , m . therefore 
the c.m. of the free neutrons has to be necessarily larger than Vk m . , that means, on av¬ 
erage, more neutron emission from the PLF side. The case of 7 = 0.8 produce instead for 
Z b > Z PL p a larger negative value (D c z ). In this case {V r z m (i ) > V c . m . and therefore the cloud 


10 











FIG. 5. (Color) For the investigated system 48 Ca + 2 ' Al we plot for different Zb windows (the 
related extremes are indicated by the double labeling in the Zf, axis) and for TKEL<350 MeV the 
measured average dipolar signal (D c z ) (red point). The experimental value are compared in the 
same figure with the results of COMD-II +GEMINI calculations (black point) filtered through the 
simulated response of the experimental apparatus and data-analysis selections. Panel A),B),C) 
show the comparisons for different 7 values characterizing the iso-vectorial interaction. The y 2 
values with respect to the experimental data are also shown. Panel D): the calculated values 
of (Dz) for 7 = 1 are compared to the experimental values for ( D c z ). The error bars represent 
uncertainties due to the statistics of the simulations and measurement. The estimated uncertainties 
related to experimental values (D z ) are in many cases smaller than the plotted symbols. 

of neutrons has a c.m. velocity smaller than V c . m .. I 11 this case we have on average more 
neutron emission from the mid-rapidity and TLF side. A closer look to the fragments put 
in evidence in this case a dominance of processes producing a PLF fragment, only partially 
equilibrated in charge/mass and a almost complete disassembly of the TLF. Up to now we 
have discussed the behavior of (D z ) which shows a good sensitivity to the parameters of 
the effective interaction. To recover information on the global degree of isospin equilibration 
we can evaluate (Dz) through calculations by using the same set of parameters which ’’best 
fit” the experimental value of (D c z ). In Fig.5 (panel D)), we compare the calculated values 
of (Dz) for 7 = 1 with CoMD-II+GEMINI (including the efficiency effect) with the exper- 
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imental values ( D c z ). The connection between the two quantities can be approximated by 
the following simple relation: (Dz) = (D c z ) + z it((V c Z mJ -Vcm) where (V c z ) is the c.m. 
velocity for the subsystem of the charged particles. Therefore the second term in the above 


expression give us an estimation of the effect associated to the undetected free neutrons 


34] 


which clearly participate in to determining the global degree of isospin equilibration. Before 
to conclude this section we add some observations about the ’’invariance” of the discussed 
quantities with respect statistical decay modes. As previously observed this was shown in a 


rather general way in ref. [7]. In the present work we have also checked if this insensitivity 
is still maintained by operating with a realistic efficiency of the CHIMERA apparatus. In 
particular, in the ’’simulated” analysis we can keep the memory of the selected primary 
events from the dynamical model ordered according to the different Z b values and TKEL. It 
is therefore possible to evaluate the same quantity from the dynamical model (D§) without 
taking into account the GEMINI secondary decay processes and comparing it to the value 
of (Dz) obtained from the complete calculations including the efficiency effects. The calcu¬ 
lations have been performed for different Zb windows and for TKEL < 350 MeV. Table I, 
as an example, collects the obtained results for different Zb bins by using the calculations 
for 7 = 1. From the table we can appreciate the close correspondence between the results 


TABLE I. For different windows of Zb and for TKEL< 350 MeV the values (cm/nsec) of (Dz) and 
the corresponding values of (D z ) are shown for the 48 Ca + 27 Al system for 7 = I 


Z b 

(D z ) (cm/nsec) (D§ 

) (cm/nsec) 

12-15 

-5.73 

-5.9 

15-17 

-8.36 

-8.34 

17-19 

- 8.86 

-8.87 

20 

-9.79 

-9.78 

21-22 

-5.36 

-5.30 


after the GEMINI 


32] de-excitation stage and the ones related to the dynamical primary 


events. 
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IV. SUMMARY AND OUTLOOKS 


In summary, dipolar degrees of freedom for the collision l8 Ca + 27 Al at 40 MeV/nucleon 
have been investigated for the first time by means of 47r multi-detector CHIMERA. The 
study has been carried out through the measurement of the observable (D z ) (associated to 
all the emitted charged particles) and through a related estimation of the (Dz) values which 
are closely linked with the global Isospin equilibration process along the beam direction. The 
discussed insensitivity to the statistical decay make this observable able to investigate in an 
exclusive manner the overall dynamics on this phenomenon. The study performed on the 
behavior of ( D c z ) and (Dz) as function of the size of the biggest fragments for TKEL< 350 
MeV shows a noticeable sensitivity of this observable to the parameters of the microscopic 
effective interaction. First attempts to reproduce the experimental values of (D z ) in essen¬ 
tially binary processes show an overall agrement with CoMD-II calculations using a stiffness 
parameter for the iso-vectorial interaction 7 ~ 1 -r 1.2. More experimental investigation 
should involve different systems and different reaction mechanisms. In particular, a next 
step forward in this kind of measurements would require a more detailed investigation includ¬ 
ing a reliable valuation and/or minimizations of possible systematic errors on the velocity 
of the produced charged particles. This would permit a reliable experimental estimation of 
(Dz) allowing also for a corresponding experimental valuation of the global effect associ¬ 
ated to the dynamically emitted neutrons. For this purpose long measurements involving 
targets and projectiles having the same charge/mass asymmetry (vanishing values of (Dz) 
independently from the reaction mechanism) are necessary. A more detailed analysis could 
also allows to estimate the specific contribution to the equilibration process produced by the 
prompt and mid-rapidity emission as compared to the transfer of charge/mass between the 
main partners. The more challenging attempt to investigate on the other component along 
the X impact parameter direction could be of relevant interest. ( Dx) in fact is linked, in a 
global way, and in an independent way from statistical emission, with the differential flow 
of particles characterized by different charge/mass asymmetries values 35]. 
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